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Abstract 

One possible evolutionary scenario of the dense gluon system produced in an ultrarelativistic 
heavy ion collision is the bottom-up thermalization scenario, which describes the dynamics 
of the system shortly after the collision via the decay of originally produced hard gluons 
to soft ones through QCD branching processes. The soft gluons form a thermal bath that 
subsequently reaches thermalization and/or equilibration. There is a scaling solution to the 
bottom-up problem that interpolates between its early stage, which has a highly anisotropic 
gluon distribution, and its hnal stage of equilibration which occurs later. Such a solution 
depends on a single parameter, the so called momentum asymmetry parameter 5. With 
this scaling solution, the bottom-up scenario gets modified and the evolving parton system, 
referred to as the m’bottom-up parton system throughout this paper, is described by this 
modification. The time evolution of the system in the original bottom-up ansatz is driven 
by the saturation scale, Qs- However, for the m’bottom-up we generalize the ansatz of 
the evolution by introducing two additional momentum scales, which give a thermalization 
time and temperature of the soft gluon bath somewhat different from those obtained when 
the m’bottom-up matches onto the final stage of the original bottom-up scenario. 
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1. Introduction 

The original bottom-up thermalization ansatz [1] in its basis has the parton satura¬ 
tion mechanism expected to be valid at high parton densities. It leads to thermaliza- 
tion/equilibration of a parton system produced after a heavy ion collision. This scenario 
emphasizes the importance of branching processes of initially produced gluons, and as a 
consequence of such processes the total number of gluons increases between the initial and 
thermalization times of the evolution of the system. 

The bottom-up does a good job in predictions of charged hadron multiplicities observed 
at RHIC and LHC energies [2]. However, it should be noted that the overall bottom- 




up picture can also be modified [3, 4] to include other possibilities for evolution of the 
produced matter, which we shall refer to as modified bottom-up parton system or simply 
as m’bottom-up parton system. The modification is done by a scaling solution to pre¬ 
equilibrium evolution which interpolates between initial plasma instabilities [5, 6, 7] present 
in the dense gluon system produced immediately after a collision, and the final equilibration. 
Depending on a single parameter known as the momentum asymmetry parameter 5, the 
scaling solution matches onto the original bottom-up picture either at an intermediate stage 
or toward the end of the evolution given by the original bottom-up. 

Thus, by having the original picture modified we further argue that in the m’bottom-up 
scenario the dominant qualitative and semi-quantitative features of the evolution of a parton 
system can be described by two momentum scales, D and Note that in the original 
bottom-up the evolution is described only by one scale - the saturation momentum Qg. 
The situation somewhat resembles the overpopulated Glasma, where two dynamical scales 
A and are introduced to describe the dynamics of the system from the overpopulated 
initial stage all the way to thermalization [8]. In this paper we show that the scales D and 
can describe the dynamics of the m’bottom-up parton system from its initial stage toward 
thermalization. However, depending on the parameter 6, the thermalization time and the 
thermalization temperature obtained, can be of the same order (or different) as compared 
with those derived in the original bottom-up scenario. The results of this paper can be 
useful for calculating the yield of direct photons and di-electrons as well as the elliptic flow 
of direct photons, in particular, in Au-|-Au collisions at = 200 GeV collision energy, 

measured by the PHENIX experimental collaboration [9, 10, 11, 12]. 

In the next section we give a brief overview of the Glasma dynamical scales. In the 
third section we derive some quantities, which describe quark-gluon matter, by using the 
dynamical scales H and of the m’bottom-up scenario. Besides, we also show that these 
and Glasma scales are the same at least on the parametric level. In the last section we 
show some estimates of the thermalization time and temperature at 200 GeV. 

2. The two momentum scales of the Glasma 

First, let us discuss what happens in the Glasma. The Glasma is theorized to exist 
from the earliest time after a heavy ion collision when the helds are considered to be highly 
coherent, and when most of the energy is in coherent field degrees of freedom rather than in 
the degrees of freedom of incoherent quarks and gluons. Thereby, the matter in this phase is 
neither the Quark-Gluon Plasma nor the Golor Glass Condensate, however, it has features 
of both. Almost instantaneously after the collision, the transverse color helds of the Color 
Glass Condensate transform into the longitudinal color electric and color magnetic helds of 
the Glasma [13]. In this scenario the gluons, in the sense of particles, are produced from 
the classical evolution of such color electric and magnetic helds, and after initial plasma 
instabilities produce a distribution that is approximately isotropic in momentum space. 
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The evolutionary and thermalization processes of the Glasma and m’bottom-up sce¬ 
narios are different from each other because of the underlying mechanisms by which they 
work. Nonetheless, they can have some properties which are similar phenomenologically, 
and we shall investigate these properties in this paper. In this section we follow the lines 
of reasoning represented in Refs. [8] and [14]. 

It is expected that on the transverse momentum and mass scales the effects stemming 
from evolution to the thermalized distribution of gluons are enhanced. The evolution of the 
system toward thermalization can be traced based on the Boltzmann transport equation: 

dtfik,x) = Ck[f], (1) 

where the Ck[f] is the collision integral, and the f{k, x) is the particle distribution function. 
The main qualitative features of the solution of the Boltzmann equation can be described 
if it is assumed that the evolution is dominated by only two scales - the “infrared” and 
“ultraviolet” A, characterizing the gluon distributions in the Glasma. It is also assumed 
that the elastic collisions play a dominant role in driving the initial gluon distribution 
toward local equilibrium. 

The scale A^ is a momentum scale at which the distributions are highly coherent, and 
it is time dependent. The scale A above which the distributions become dilute is also time 
dependent, which coincides at the earliest time with A^: A(to) = As(ro) ~ Qs (the r is 
the proper time of the collision). Along with time these scales separate from each other 
where the A^ decreases rapidly, and the A evolves more slowly. Upon reaching equilibration, 
the A becomes the initial temperature of the Quark-Gluon Plasma: A.{Ttherm) ~ Tin^QCP- 
On the other hand, the infrared scale becomes the so called non-perturbative “magnetic 
scale” [15, 16] in the thermalized matter: ^s{Ptherm) ~ OisTth^rm = otsTin^gcp-, where the 
as is the QCD coupling constant. The thermalization is accomplished by parametrically 
splitting apart these initially overlapping momentum scales by and the corresponding 
time is determined from the following requirement: 


^s{Ptherm) ~ CXs-^iPtherm) ; ( 2 ) 

at which the gluon occupation number, /(A), becomes of the order of unity. Generally, 
the /(A) is proportional to ~ l/a^ which is down to the ultraviolet scale, whereas the 
thermal distribution, l/cop, builds up from the infrared scale gradually. Thus, the gluon 
distribution becomes a thermal gluon distribution function when the infrared scale satisfies 
the condition in Eq. (2). 

In the momentum space there are the following cases (at time r > 1/Qs) in the corre- 
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spending three ranges: 


fgip)^— atp<As, 

as 

fg{p) ^ at As <p < A, 

as Up 

/g(p)~0 atp>A, (3) 

where the p is the gluon momentum, and the Up its energy. In general, the occupation 
number is expressed as 

= ( 4 ) 

Additionally, the gluon density and the Debye mass are expressed via the scales A and A^; 

Ng ~ —A^A,, (5) 

as 

Ml ~ AA,. (6) 

At the initial time, by having the A(ro) and As(ro), these two equations are represented as 


AT 


Ml^Ql ( 8 ) 

At the thermalization time they will have forms represented by the initial temperature of 
the Quark-Gluon Plasma: 

^g ~ '^in,QGP ) ( 9 ) 

Ml ~ oisTI qqp . ( 10 ) 

The time evolution is dominated by the gluon density, and there can be some fixed asym¬ 
metry between the typical transverse and longitudinal momentum scales characterized by 
a parameter 5' which is defined in terms of the longitudinal pressure: 

PL = 5'e, ( 11 ) 

where 0 < 8' < 1/3, with d' = 0 corresponding to the maximal anisotropy between 
the longitudinal and transverse pressure, and 8' = 1/3 corresponding to the isotropic 
expansion. Finally, since the electromagnetic particle production arises from quark charges, 
the quark distribution function should also be included in the overall evolutionary picture: 

fg=Fgip/A), (12) 
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that brings out the quark density as 


Nq^A^. (13) 

Initially not many quarks are present in the matter and this is suppressed by the power 
of as, which is small if the Qs is large compared to the QCD scale Aqcd- Namely, at 
the earliest times, Nq ~ agNg <C Ng, nevertheless, later the two densities approach each 
other becoming parametrically of the same order, Nq ~ Ng. Therefore, in the time interval 
1/Qs ^ T Ttherm the quark density increases to a value being of the order of the gluon 
density, and is no longer suppressed. Thus, the quark production at r ~ l/Q^ might not 
be significant, or anyway not as significant as at later times. 

3. Prom the original bottom-up thermalization to the m’bottom-up thermaliza- 
tion 

3.1. The scaling solution and the matching onto the original bottom-up 

Now let us discuss what happens in the m’bottom-up scenario. The bottom-up evolution 
is also based on using the solution of the Boltzmann equation which, with inclusion of the 
particle production, makes the parton system approach the kinetic equilibration during 
a time of the order of r ~ as obtained in [1]. The thermalization occurs in 

the limit of Qs Aqcd , which corresponds to very large nuclei or very high collision 
energy. Let us start with the original bottom-up, which is a mechanism that describes 
the system evolving toward the thermalization/equilibration after passing through three 
distinct stages. 

Stage 1. Shortly after the collision, hard gluons are produced being distributed highly 
anisotropically. These gluons dominate in the time range of 1/Qs < t < as 

Stage 2. In this phase these initial hard gluons still dominate but start emitting soft 
gluons via QCD branching. The soft gluons in turn equilibrate and form a thermal 
bath, which initially carries only a small fraction of the total energy of the system. 

Also, the thermal bath draws energy from the hard gluons. These processes occur 
— 3/2 — 5/2 

in the range of a* ' < QsT < as . The soft gluons start to overwhelm the hard 

gluons, in terms of number, at t ^ as 

Stages. Now the soft gluons are thermalized but continue drawing energy from the 
hard gluons until all the energy is removed from them and the whole system is equi¬ 
librated. The full thermalization is achieved when the primary gluons have lost all 
their energy to the soft gluon bath. Parametrically, this occurs at r ~ as 
which is the end of the time range of this final stage: as < QsT < cts . 
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However, after taking into account the existence of collective effects in the form of 
plasma instabilities [5, 6, 7] in the initial stage of the bottom-up picture, one can think 
about a rapid thermalization scenario. It is because the plasma instabilities occurring 
in the dense gluon system produced immediately after a collision, increase the amount 
of energy transformation from initially produced hard modes to soft modes radiated by 
the hard ones. Meanwhile, it has been shown in [17] that the full equilibration time in 
the presence of the instabilities is not much shorter relative to that of the bottom-up. 
But these instabilities cannot lead directly to the equilibration since they would give an 
equilibration time parametrically of the order of QsT ~ 1. For a more detailed insight into 
this problem, a scaling solution has been proposed [3, 4] for following the evolution between 
the instabilities taking place in the initial phase and the system in the final equilibration. 
This solution, depending on one single parameter 6, matches onto the intermediate stage 
and/or the late stage of the evolution of the system given by the bottom-up thermalization, 
where the 5 accepts values in the range of 0 < <5 < 10/21. As in the case of the Glasma, 
here the 5 again describes the asymmetry between the transverse and longitudinal scales of 
the gluon interactions. The scaling solution makes the bottom-up scenario getting modified 
to the more general m’bottom-up scenario. The number 10/21 is the absolute limit the 5 
may accept in the m’bottom-up. So the proposed scaling solution is given by the following 
set of equations: 




O^sfs 


Ql 


Qs 


asiQsT)^^^ ^ ’ 

1 

" (Q.t)V3+5/5 ’ 


Md 


(Q,r)V3-25/5 ’ 

Qs 

(Q.r)V2-35/10 


(14) 


where the Ng is the number density of the soft gluons, the kg is the soft gluon momentum, 
the fs is the soft gluon occupation number, and the Md is the Debye mass. 


3.2. The two momentum scales of the m’bottom-up scenario 

There might be a question about whether or not similar (or different) momentum scales, 
such as the discussed Glasma A and A^ scales, are applicable for the m’bottom-up. Con¬ 
ceptually it should take place, since as it was stated before the main qualitative features 
of the solution of the Boltzmann equation can be described with an assumption that the 
evolution is dominated by such scales. In the case of the m’bottom-up let us designate 
these scales as D and and check whether or not such scales can really exist. Making use 
of the second formula of Eq. (3), we assume that an equivalent relation (with p ~ Wp at 
Mo <C ujp) can exist in the m’bottom-up picture, though it is not the case in the original 
bottom-up: 

f(p) ~- - at ^Ig < p < Q . 

ag p 
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(15) 







These scales are related to the soft gluon density Ng-. 


Ns= ra^pfip)-^ rdp—— 


Us p as 


and to the Debye mass 


rO. 


nf2 I ^3 fiP) /"^ 2 j ^ on 

Mjy = Us d p - ^ as p dp -^ ~ UUs 

J P J asp^ 

For the quark occupation number we use /g ~ 1, such that 

Nq= r d^pfip)r^n^. 


(16) 


(17) 


(18) 


As in the case of Glasma, this also comes from the reasoning that at the earliest times 
we have Nq ~ agN^ ^ Nh but at later times (for example, at thermalization) Nq ^ Ng. 
Here, the is the number density of the primary hard gluons (initial gluons). 

So how can we make sure that the scales D and really hold in the m’bottom-up 
scenario ? One may see this if we look back at Eq. (5) and Eq. (6), and start with the 
analogical equations for the m’bottom-up such as obtained in Eq. (16) and Eq. (17): 

Ng ~ —n^ng , (19) 

ag 

Ml ~ nVlg , ( 20 ) 

with the following requirements: 

D(to) = Hs(ro) ~ Qg , (21) 


and 


^{j~therrri) ^therm ^ ; ^sirtherm) ~ ®s^(A/ierm) ■ (^^) 

Thereby, using Eqs. (19), (20), (21) and (22) at the initial time tq ~ Q^^-, we shall have 


the “initial conditions” 


Nh 


mI ~ 


ag{QgTQ) 

Ql 


(23) 


(24) 


{QsTq) ’ 

as it holds in the original bottom-up ansatz as well. At the thermalization time there should 
be the following conditions, which we name as “boundary conditions”: 


~ , (25) 

Ml ~ agTl^^j.^ . (26) 
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3.3. The derivation of the soft gluon number density, Debye mass, energy density and 
entropy density via thermalization temperature 
In this section we show how Eqs. (25) and (26) can be derived in the m’bottom-up. For 
this purpose we should use the Ng and Mo from Eq. (14), along with the thermal bath 
temperature T{t) from [3, 4]: 


Ng 


asiQsT)^^^ ^ ’ 


M] 


Q 


D 


(Q^^)1-35/5 ’ 


~ q2 ^2(35-78(5)/(39(5-10)|'q ^^2(15-36(5)/(39(5-10) 


(27) 


(28) 

(29) 


Independent of what value the 5 takes, as long as S > 1/3, the temperature in Eq. (29) 
and the scaling solution in Eq. (14) match up onto the final stage of the original bottom-up 

only at the final time QgT ~ .In this case the temperature in Eq. (29) reduces to 

2/5 

T ~ QsOis which is independent of <5. 

Consequently, we should search for some solutions for Ng and Mo that should be 6- 
independent. We start with Eq. (28) and Eq. (20) such as 


Ml ~ 


Ql 


and solve them for 


{QgTy-35/5 


Tig 


and 


Ml ~ 


TlTlg 


Q 


(30) 


(31) 


D (Q,r)i-3^/5 • 

But still the scale 12 (dependent on r) stands in the denominator of Eq. (31). Then making 
use of Eqs. (27) and (19) 


Ng 


Q 


and Ng - D^Dg, 

(Xg 


as(Qs7-)‘^/3 <5 

along with (31) and solving them for Tl, we further find the scale as a function of r: 

^ ^ (5-6<5)/15 


(32) 


n ~ Qo 


QsT 


(33) 


Thereby, by inserting the Tl from Eq. (33) into the denominator of Eq. (31), we obtain the 
scale Tig as a function of r. 

^ ^ (10-35)/15 


QsT 


Tig ~ Qg 


(34) 









~n, ~n, 



Figure 1: The two scales Q and fls as functions of the proper time t, at three values of 5 shown from top 
to bottom as (5 = 0, h = 1/3 and S = 10/21. 


In Fig. 1, the derived time-dependent scales 1/ and are parametrically depicted at selected 
three values of the 6 parameter. 

With these two derivations of and we can parametrically find the thermalization 
time in the m’bottom-up scenario. We can get it using the condition in Eq. (22): 


Qs 


1 


Q 


sTtherm 


(10-3(S)/15 

~ asQs 


1 


Q 


sTtherm 


(5-6(5)/15 


(35) 


which ultimately gives 

QsTtherm ~ . (36) 

As a next step we divide both sides of Eq. (27) and Eq. (29) (with T^) on each other, 
and both sides of Eq. (28) and Eq. (29) again on each other (also by using QsTtherm from 
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Eq. (36)). Then we shall have the following two formulas: 


Ns 


r\j 


r\j 




[1+3(35-785)/(395-10)] / -15/(5+35)^ 

I therm 

[-1955+50-11752+305-525-315<5+1170(5+70252]/[(5+35)(395_io)] 

CKg A 

[675-1620<5+780(S-200-58552+1505]/[(5+3(5)(39<5-10)]^3 
^ therm 

[585<52-585(52+690<5-690(5-475+475]/[(5+3(5)(39<5-10)]^3 

therm ’ 


( 37 ) 


which reduces to Ng ~ 7)^erm> same as Eq. (25). And 

2 ^ „-2(35-785)/(39^-10) / -15/(5+35)^ 

iViD CXg I I therm 

^ [-350-210(5+780(5+46852+585(5-150-35152+90<5+450-10805]/[(5+3<5)(395-10)]^2 

therm 

^ [11752+l655-50]/[(5+35)(395-10)]Tn2 

^ therm ’ 


(38) 


which reduces to M|, ~ the same as Eq. (26). Alternatively, an analogous 

outcome can be demonstrated if we use the O and scales. 


A^ 

r3 


\therm 


Ns 


and 


Ml, 

2^2 


\therm 

Ml 


/ (q 3 ^3(35-785)7(395-10) (g^^)3(15-365)/(395-10)^ 


2^ 

^ therm ’ 


(39) 


(Sia.) / [qI 


OsT. 


therm ’ 


(40) 


Consequently, the results in Eqs. (37) and (38) are independent of any value the 5 accepts. 
Note that the thermalization temperature Tt^erm of th® m’bottom-up parton matter is 
assumed to be parametrically of the same order as the initial temperature of the Quark- 
Gluon Plasma Tm.QGP- 

However, this is not the end of the story. An equilibrated system of soft gluons at the 
thermalization temperature Ttherm must satisfy the condition ~ for the energy 

density. Needless to say that we should also be able to prove the validity of this relation 
in our ansatz. From [3, 4] we hnd the energy density of the soft sector expressed via the 
soft gluon number density and momentum, ~ Ngkg. Then from the scaling solution of 
Eq. (14) we will have the following expression: 




a.(Q.r)(25-215)/15 ’ 


(41) 
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By dividing both sides of Eq. (41) and Eq. (29) (with T^) on each other along with using 
QsTtherm from Eq. (36), the energy density becomes 


~ a 


- [l+4(35-785)/(39(5-10)] (^-15/(5+3(5) 


-[(25-215)/15+4(15-36(5)/(395-10)] 


rpA 

therm 


~ a 


-195<5+50-117<52+30<5-700-420(5+1560<5+936(52]/[(5+3(5)(39<5-10)] 


~ Oi 


[975<5-250-81952+210<5+900-2160(5]/[(5+35)(39(5-10)]Tn4 

A Cig -L ih 

[0+0+0]/[(5+35) (395-10)] t^4 

c +h 


therm 


therm ' 


(42) 


In addition to the above derivations, one can also pay attention to the following ob¬ 
servation. Let us assume for a moment that the scale fl from Eq. (33) and the scale 
A ~ Qs {1 /Q from [8] are parametrically the same. 


n ~ A 


Q. 


QsT 


(5-65)/15 


Qs 


QsT 


(l+25')/7 


Mathematically this assumption is valid if 


or vice versa 


5 = 


5' = 


10 - 155' 
21 

10-21(5 

15 


(43) 


(44) 


(45) 


Then what will happen if this (5 ++ A relation is being inserted into the r.h.s. of Eq. (34) ? 
This results in what follows: 


^ Qs 


1 


QsT 


(10-35)/15 


Qs 


1 


QsT 


(4+5')/7 


(46) 


But the r.h.s of Eq. (46) is actually the scale A^ from [8], which is A^ ~ Qg (l/Q^r)^^"*"'^ 

In this connection let us again take a look at Ref. [8] but for the gluon density and Debye 
mass 


AT,- 


Q 


Mi 


Q 


a.(g.r)(6+5<5')/7 ^ (Q^^)(5+35')/7 ’ 

as well as for the thermalization time, which is obtained with the condition of Eq. (2) 

QsTtherm ~ Ot~'^^■ 

Besides, we also wish to take a look at the energy density in gluon modes 

Eg ~ —AsA^ , 

Ois 


(47) 


(48) 


( 49 ) 
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which, in the case of the longitudinal expansion under the assumption of the parameter 5' 
independent of time, is represented as an equation for the evolution of the energy density: 

~ eg(To) j > (50) 

Thereby, by making use of Eq. (45) it is easy to check that the following transformations 
take place: 


Ng of Eq. (47) 

—)> 

Ng of Eq. (14) ; 

(51) 

Md of Eq. (47) 


Md of Eq. (14) ; 

(52) 

rtherm of Eq. (48) 

—)• 

'^therm of Eq. (36) , 

(53) 

eg of Eq. (49) 


e, of Eq. (41) . 

(54) 


However, these transformations can be valid in the range of 10/21 > 6 > 5/21 correspond¬ 
ing to the range of 0 < <5^ < 1/3, This can be seen as follows: 

at (5' = 0 ^ = 10/21, at 5' = 1/3 ^ = 5/21, at 5' = 2/3 ^ 5 = 0. 

On the other hand, if formally one needs to recover the static case by setting 6' = — 1, 
corresponding to constant energy density, eg(r) —eg(ro), then it will also be the case for 
the m’bottom-up scenario, €s{t) —)■ €s{to) ~ Q^/os, because in that case 5 = 25/21. 

For the energy density one may simultaneously see 


Eg ~ — 
CXa 


Qt 


a.(Q.r)(25-2W)/i5 


(55) 


As a last step we wish to prove that the condition for the entropy density in the thermal 
gluon bath at thermalization is derived in the form of s ~ '^fherm ■ 


T3 


\therm 


HV ('Q3^3(35-785)/(395-10)(g^.^)3(15-365)/(395-10)^ 


a„ 


-3(35-785)/(395-10) /^-15/(5+3^)r 


therm 


[-525-315(5+11705-I-702<52+5855-150-702<52+1805-I-675-16205]/[(5+35)(39<5-10)]t-.3 

o +h 


^ Jo+0+0]/[(5+35)(39(5-10)]j.3 

therm 


therm 


(56) 


And 

solution 


at the end of this section we mention one of the results from [3]. If 
in Eq. (14) changes the character at a time r given by 

Q.r ~ , 


5 > 1/3, the 


(57) 
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at which /^ ~ 1. In that case the scaling solution goes into evolution much like the final 
phase of the original bottom-up where the soft gluons are thermalized and the harder 
gluons feed energy into the soft thermalized system causing the temperature to rise with 
the time until the whole system is thermalized. We see that Eq. (36), which we get after 
introducing the scales O and into the m’bottom-up scenario, turns out to be the same as 
Eq. (57), though in our approach the corresponding time is already considered to be as the 
thermalization time of the system stemming from the thermalization condition of Eq. (22). 


4. The thermalization time and the thermalization temperature of the m’bottom- 
up parton system at RHIC = 200 GeV 

As an example of our formalism, one may find values of the thermalization time (from 
Eq. (36)) and temperature (from Eq. (29)) at RHIC = 200 GeE collision energy by 

following the procedure that has been used in [2]. So for the thermalization time we have 

Ttherm = Ctherra Q~^ , (58) 


where Ctherm is the so called thermalization constant. Eor the thermalization temperature 
we have 

= Ct . (59) 

which by using Eq. (58) reduces to 

Ttherm ^ 0.16543 C a(5-65)/(5+35)g^ ^ 

where we also use the numerical constant Ct expressed by the “gluon liberation” coefficient: 

(61) 


0.16543 C. 


The C is a parameter that links the number of gluons in the nucleus wave function to the 
number of gluons which are freed during the collision. 

As discussed already, the number of gluons increases with the time r, because the hard 
gluons degrade and the soft ones are formed and start to dominate in the system. Such an 
increase of the number of gluons must be 2 or larger, and can be found as the following 
ratio. 


R = 


[A^s(t)((5sT)] 


'^the 


> 2 


[A^/x('r)((55T)] 1^0 

R ~ 0.13061 > 2 


2 < 0.13061 I -- 
9 ) 




(15(5-5)/(5-l-35) 


A2 

/'■QCD^ 


(62) 
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where we make use of Eq. (58) and Eq. (60) along with the formulas shown in what follows: 


Ns{Ttherm) — 2(A^c 1) 2 '^therm 1 

'TT'^ 


TT^ 


A'h(ro) = C 
as{Ql) - - 


dvr^iVcas {QsTq) ’ 
dvr 


(11 

For Nc = 3 we take Nf = 3. Also, by using the well-known relation at midrapidity 

A/(l-|-A/2) 


Qs{A, \/i) — Qq{A, -v/io) 


KV^oJ 


(63) 

(64) 

(65) 


( 66 ) 


with A = 0.288 [18], we shall have the saturation momentum ~ 1.115 at ^/s]^ = 
200 GeE in Au-|-Au collisions, which is obtained from Qq = 1 GeV^ at -s/snn = 130 GeE 
that has been used in [2]. Thus, the ratio R in Eq. (62) can be 2 or larger if the overall 
constant G^ taken adequately. 

We need also one more formula with G and Ctherm after which these parameters can be 
determined. That formula is derived by comparing the calculable charged hadron multiplic¬ 
ity at midrapidity at ^Jsnn = 200 GeE with that from RHIC Au-|-Au data. As a reference 
value we use the result by the PHOBOS collaboration [19], namely 3.78±0.25 (syst). Thus, 


2 dNch 


'N, 


part 


dr] 


)\exp = 3.78 ± 0.25 {syst). 


The charged hadron multiplicity can be calculated as 


2 dNch I R C 


N, 


part 


dr] 


In 


Q 


A2 

^'■QCD ^ 


which is further simplified to be 
2 dNch \ 


Npart dr] 


'dvrA E-155)/(5+3<5) 
~ 0.04354 ( j 


^3 ^(35-695)/(395-10) 
U O,; 


''therm 


In 


Ql 


(67) 


( 68 ) 


18<5/(5-I-3(5) 


A2 

^'-QCD 


, (69) 


where we use Eq. (62) for the ratio R. Here the charged hadron multiplicity will be equal 
to 3.78 if the overall constant G^ taken adequately. 
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Let us now find the Ttherm and Ttherm as functions of the 5 parameter. We will try two 
values for the gluon liberation coefficient: the first one calculated in [20] with C = 21n2 ~ 
1.386, and the second one in [21] with C = 1.1. For example, at <5 = 1/3 we will have 


Ctherm — 


2.260 

C3/4 


and Ctherm ^ 


1.978 

C2/4 


^ C < 1.704 and Ctherm > 1-515. (70) 

If we use the above values of the “gluon liberation” coefficient, then from Eq. (58) and 
Eq. (60) we shall have the following fixed values for the thermalization time and tempera¬ 
ture: 


At (7 = 1.1 Ttherm - 3.45 fui , Ttherm - 2Q2MeV. 

At C = 2 In 2 Ttherm - 2.90 fm , Ttherm - 277 MeV. (71) 


Also, at C = 1.1 =► i? ~ 3.1, and at (7 = 2 In 2 =► ii~ 2.46. In Eig. 2 and Fig. 3 
the solid lines show the d-dependent thermalization time Ttherm at both values of C. The 
(5-dependent thermalization temperature Ttherm at both values of C is shown in Fig. 4. 

One may put some constraints on the possible values of the parameter 5. It can be 
done if we find a time ri at which Nh = Ng- Namely, by using Eq. (64) and Eq. (63) with 
Eq. (60) we derive 


QsTl 


( - - - 

V0.1306 (72 


(39(5-10)/(35-69(5) 

-(95-195(5)/(35-69(5) 

(Jig 


(72) 


It is obvious that at 5 = 1/3, this formula becomes parametrically the same as that of the 
original bottom-up, i.e., Nh = Ng at QgTi ~ Os ' . Then if we calculate the time ri at 
C = 1.1 and (7 = 2 In 2 for 5 = 1/3, we shall obtain the following result: 


At (7 = 1.1 —> Ti~2.60/m, 

At (7 = 2 In 2 ^ ri~2.31/m. (73) 


In Fig. 2 and Fig. 3 the dashed lines show the 5-dependent ti at both values of (7. So 
since the soft gluons start to overwhelm, in terms of number, the primary hard gluons at 
Ti ~ Ug i95'^)/(35 69(5) means that the following inequality must always take place: 
Ttherm > Ti- By Comparing the results shown in Fig. 2 and Fig. 3, we see that the cases of 
5 < 0.26 should be excluded. On the other hand, we notice that the evolutionary picture 
looks highly unlikely at larger values of 5 close to 5 = 10/21, which shows that the system 
never gets thermalized. Thus, we find that the thermalization time and temperature at 
values of 5 from a range around 5 = 1/3 is comparable to those from [2]. 
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At C=l.l 


^therm (fttl) 



Figure 2: The solid and dashed curves show the (5-dependent Ttherm and ti at C = 1.1, respectively. 


At C=21n2 

ftherm (fttl) 



Figure 3: The solid and dashed curves show the 5-dependent Ttherm and ri at C = 2 In 2, respectively. 


However, we should note that the results in Figs. 2, 3 and 4 are approximate because 
we make use of the parameter Ct from Eq. (61) (by having it from [1, 2]) in our derivations 
of this section. It is possible that the Ct can be d-dependent, which perhaps will (or will 
not) make the curves less steeper at large values of 5, and/or can alter the results for Ttherm 
and Ttherm to be similar to those from Ref. [2] at small values of 8. But the derivation of 
the 5-dependent Ct parameter is beyond the scope of this paper, which however will be 
considered in another work. 
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The dashed line is at C=21n2 


The solid line is at C=l.l 


rtherm (GeV) 



Figure 4: The solid and dashed curves show the (5-dependent Ttherm at (7 = 1.1 and (7 = 2 In 2, respectively. 
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